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ABSTRACT 

We present the discovery of a new quadruply lensed quasar. The lens system, SDSS 
J1330+1810 at Zs = 1.393, was identified as a lens candidate from the spectroscopic 
sample of the Sloan Digital Sky Survey. Optical and near-infrared images clearly show 
four quasar images with a maximum image separation of 1'.'76, as well as a bright 
lensing galaxy. We measure a redshift of the lensing galaxy of zi = 0.373 from ab- 
sorption features in the spectrum. We find a foreground group of galaxies at z = 0.31 
centred ~ 120" southwest of the lens system. Simple mass models fit the data quite 
well, including the flux ratios between images, although the lens galaxy appears to be 
~ 1 mag brighter than expected by the Faber-Jackson relation. Our mass modelling 
suggests that shear from nearby structure is affecting the lens potential. 
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(SDSS J133018.65-hl81032.1) 



1 INTRODUCTION 

Thus far about 100 gravitationally lensed quasars are 
known, of which ~ 30 are quadruple (four-image) lenses. 
The number ratio of quadruple lenses to double (two-image) 
lenses contains information on both the shapes of lens- 
ing galaxies and the lum i nosity funct i on of source quasars 
(Rusin & Togmark' 2001'; 'Chae' '2003"; 'Hutere r et all l2005l : 
[Oguri2007bi : ,Mandclbaum et al, 200 8). In addition, quadru- 
ple lenses allow more detailed mass modelling of indi- 
vidual lenses. For instance, the larger number of images 
provides more constraints on the lens potential, which is 
essential in probing the effects of external perturbations 



* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. Based 
on observations obtained with the Apache Point Observatory 3.5- 
meter telescope, which is owned and operated by the Astrophys- 
ical Research Consortium. Use of the UH 2.2-m telescope for the 
observations is supported by NAOJ. 
f E-mail: oguri@slac.stanford.edu 



on primary lenses (jKeeton et al.l Il997l ) and constraining 
the Hubble constant fro m time delay measurements (e.g., 
ISuvu fc Blandfordll2006h . Magnifications of merging image 
pairs in quadruple lenses satisfy distinct relations if the 
lens potential is smooth, but small-scale structures near 
the image can violate the relations. Thus flux ratios of 
quadruple lens images serve as uni que probes of substruc- 
ture or microlensing in lens galaxies (iMao fc SchneideJl998l: 

iTTTrrTrT-TTrT:^ Ti^TiTrn. r^rfrn Unnni. It^_,_i ». t^._i,_.--iJ 



iMetcalf & Madau' '2001'; 'Chiba' '20021: iDalal fc Kochanekl 

^02; Schcchtcr & Wambsganss 2002). 

In this paper, we present the discovery of a new 
gravitationally lensed quasar with four lensed images, 
SDSS J133018.65-hl81032.1 (SDSS J1330-hl810). It was dis- 
covered as part of the Sloan Digital Sky Survey Quasar Len s 
Search (SQLS: lOguri et al.ll200d l2008l : llnada et all [20081 ) . 
which ta kes advantage of the l arge spectroscopic quasar cat- 
alog fsee lSchneider et al.ll2007l ) of the Sloan Digital Sky Sur- 
vey (SDSS: Ivork et al.ll2000l ) to locate new lensed quasars. 
We place particular emphasis on mass modeling and inves- 
tigation of the structure around the lens. 
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Figure 1. The SDSS i-band image of the SDSS J1330+1810 field. The pixel scale of the image is 0'.'396 pixel -\ and the seemg was 
I'.'O. The inset at lower left shows an expanded view of the lens system. Several bright galaxies in the field are indicated by G1-G3. The 
i-band Petrosian magnitudes of these galaxies (without correcting for Galactic extinction) are 17.72 (Gl), 18.20 (G2), and 18.59 (G3). 
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Figure 2. The SDSS spectra of SDSS J1330+1810 (upper) and 
the nearby galaxy G2 (lower). See Figure [Tl for the relative po- 
sitions. The spectral resolution is i? ~ 2000, but the spectra are 
smoothed with a 5-pixel boxcar to suppress noise. Thick vertical 
line segments in the upper panel mark Mgll and Ca absorption 
lines, which suggest the lens redshift of zi = 0.373. A Mg/Fe ab- 
sorption system at z = 1.054 is also shown by thin lines. The 
feature at ~ 5600A in the spectrum of G2 is due to poor subtrac- 
tion of a strong sky emission fine (5575A). 



The outline of this paper is as follows. We describe the 
SDSS and follow-up data in Sections [2] and [3] respectively. 
The environment of the lens is discussed in Section |4] Sec- 
tion [5] is devoted to mass modelling. Our results are sum- 
marised in Section [6l Throughout the paper, we adopt the 
standard Lambda-dominated flat universe cosmology with 
matter density SIa/ = 0.26 and the Hubble constant h — 0.72 
(,Dunklev et al.„2008i ). 



2 THE SDSS DATA 

The lens system SDSS J 1330-I-1810 was first identified i n 
the data of the SDSS-II (|Adelman-McCarthv et"all 120081 '). 
The SDSS-II is a survey to map 10,000 square degrees 
of the nor thern sky with a dedicated wide-field 2.5-meter 
telescope (|Gunn et al.l l2006f ) at the Apache Point Obser- 
vatory in New Mexico, U SA. It consists of a photometric 
survey (Gunn et al.' 'l998j) with five broad-band optical fil- 
ters (Fukugita ot aL 19961 ) and a spectroscopic survey of 
quasars and gal axies selected b y a se ries of target selec- 
tion algorithms; [Richards et al.l (|2002|) present the S DSS 
quasar selection technique, and iBlanton et al.l l|2003l ) de- 
scribe the generation of the final SDSS spectroscopic tar- 
gets. The homogeneity and good quality of the data, with 
an as trometric accura cy better than O'.'l rms per coordi- 
nate (|Pier et al.l |2003| ) and a photometric zeropoint accu- 
racy better than 0.02 magnitude over the entire survey 
area (|Hogg et al.ll200ll : ISmith et al.ll2002l : llvezic et al.ll2004i : 
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Figure 3. Images of SDSS J1330+1810. All the panels are 7" X 7" in size. Nortli is up and East is left. The top row shows original 
follow-up images, the second row displays images after subtracting the lensing galaxy Gl, the third row contains images after subtracting 
4 PSFs (quasar images), and the bottom row shows the results after subtracting the galaxy and all four quasar components. Images for 
each band are displayed with the same gray-level stretch. 



Table 1. Relative astrometry and photometry from the follow-up images (see Figure[31l. The positive directions of x and y are North and 
West, respectively. The J2000 coordinates of component A are (RA, Dec) = (202.57784, 18.17562). The positional errors are estimated 
from the scatter between the five follow-up images. The errors on the magnitudes are statistical errors only and do not include systematic 
errors coming from uncertainties of PSFs and the galaxy profile. The magnitudes have not been corrected for Galactic extinction. The 
magnitudes of galaxy G refer to the total magnitudes. 



Name 


X 


y 


V (UH88) 


J (Magellan) 


H (ARC3.5m) 


H (Magellan) 


Ks (Magellan) 




[arcsec] 


[arcsec] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


A 


= 


= 


19.02 ±0.02 


18.36 ±0.05 


17.32 ± 0.05 


17.40 ±0.07 


17.06 ±0.04 


B 


0.42 ± 0.03 


-0.01 ±0.03 


19.72 ±0.03 


18.68 ±0.07 


17.54 ±0.06 


17.73 ±0.09 


17.30 ±0.05 


C 


1.30 ±0.03 


1.19 ±0.03 


19.89 ±0.01 


19.12 ±0.02 


18.18 ±0.02 


18.24 ±0.02 


17.90 ±0.03 


D 


-0.24 ±0.04 


1.58 ±0.04 


21.45 ±0.04 


19.83 ±0.05 


19.36 ±0.12 


19.13 ±0.05 


18.62 ±0.06 


G 


0.24 ±0.03 


0.97 ±0.03 


19.48 ± 0.04 


16.77 ±0.01 


16.00 ±0.01 


16.08 ± 0.02 


15.26 ±0.01 



iTucker et alll2006l : IPadmanabhan et al. 1 120081 ). are essential 
for various statistical studies. 

The SQLS identifies gravitationally lensed quasar can - 
didates using a well-defined algorithm ijOeuri et al.l l2006l ) 
appUed to the spectroscopic SDSS quasars. The algorithm 
has two parts: morphological selection, which identifies 
quasars that are poorly fitted by the Point Spread Function 
(PSF), and colour selection, which examines objects near 
each spectroscopic quasar and selects those with colours sim- 
ilar to the quasar as lens candidates. These two selections 
are designed to locate small- (~ 1") and large-separation 



(> 3") lensed quasars, respectively. The SQLS has already 
discovered > 20 new gravitationally lensed quasars from 
the SDSS data using this algorithm, includi ng both double 
and quadruple lenses (e.g., 'Kayo et al. !l2007l, and references 
therein). 

The gravitational lens SDSS J1330+1810 was selected 
as a lens candidate by the morphological selection algo- 
rithmQ Figure [T] shows the SDSS z-band image of the 



1 The SDSS photometric designation of SDSS J1330+1810 is 
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system (seeing of I'.'O). The enlarged image clearly indi- 
cates that the system, which is classified as a quasar at 
z = 1.393 from the SDSS spectrum (see Figure [2]), is not 
a point source but consists of multiple components. The 
morphology of the quasar in the SDSS image is similar 
to that of SDSS J0924+0219 l|lnada et al.ll2003l ). suggest- 
ing that this is likely to be a fol d-type quadruple l ens. 
We find no b right radio (FIRST; ' Becker et all liggsl ) or 
X-ray (RASS; IVoges et al.l 11999 ) source in the vicinity of 
SDSS J1330-I-1810. In the SDSS spectrum we find absorp- 
tion fines at ~ 3850A and ~ 5400A, consistent with Mgll 
and Ca lines due to a galaxy a,t z — 0.373. We interpret 
these features as absorptions by the lensing galaxy. In ad- 
dition a Mg/Fe absorption system &t z = 1.054 is seen in 
the spectrum. Since its redshift is quite close to the source 
redshift, the effect of the absorber a.t z = 1.054 on the lens 
potential is expected to be small. 

The SDSS image shows a concentration of red galax- 
ies ~ 120" southwest of SDSS J1330-hl810 (see Figure [1 
around Gl and G2). One of the gala xies, G2, was targeted 
by the luminous red galaxy program (lEisenstein et al]|200ll ) 
and has been observed spectroscopically in the SDSS (Figure 
[U, with a redshift of z = 0.3126. The presence of the [Oil] 
emission line (rest-frame equivalent width of — 9A) suggests 
ongoing star formation activity in galaxy G2. We explore 
the possible group further in Section |4] 



3 FOLLOWUP OBSERVATIONS 

3.1 Optical and Near-Infrared Imaging 

Near-infrared {J H K s-ha.nds) images were taken using 
Persson's Auxilliar y Nasmyth Infrared Camera (PANIC; 
iMartini et al] |2004| ) at the 6.5-meter Magellan I (Walter 
Baade) telescope on 2007 July 5. The total exposure time 
was 405 sec in J, 540 sec in H, and 360 sec in Ks- The 
seeing was 0'.'74 — 0'.'82. We also obtained _ff-band images 
with the Near-Infrared Camera and Fabry-Perot Spectrom- 
eter (NICFPS) at the Astrophysical Research Consortium 
3.5-meter (ARC3.5m) Telescope on 2007 June 1. The total 
exposure time was 1200 sec and the seeing was 0'.'88. Optical 
imaging was conducted with the Tektronix 2048x2048 CCD 
camera at the University of Hawaii 2.2-meter (UH88) tele- 
scope. A 500 sec image in 1/-band was taken on 2008 March 
6, under O'.'Sl seeing. All the data were reduced using stan- 
dard IRAF tasks. The zero-point magnitudes of the infrared 
images were est imated using Two M icron All Sky Survey 
(2MASS) data l|Skrutskie et all 120061 ') . whereas the UH88 
image was ca librated by the standard star PG0918-I-029 
(|Landoltlll993 ). 

The images shown in Figure [3] confirm that the sys- 
tem is indeed a typical fold-type quadruple lens, with two 
merging images and the other two images on the other side 
of th e lens. For defini teness, we fit the images using GAL- 
FIT (|Peng et al.ll200^ . The assumed model consists of four 
PSFs and a lens galaxy modelled by a Sersic profile with 
the convolution of the PSF. We adopt nearby stars as PSF 



5308/40/1/73/88 (run/rerun/camcol/filed/ID), and the spectro- 
scopic designation is 2641/255/54230 (plate/fiber/MJD). 



templates. We first left the Sersic index n as a free param- 
eter and found the best-fit value to be n ~ 3.4, which is 
close to the canonical value for early type galaxies, n = 4. 
Thus in what follows we fix the Sersic index to n = 4. We 
find that this model fits the data quite well. The subtracted 
images shown in Figure [3] show virtually no residuals. Ta- 
ble [1] summarises the relative astrometry and photometry 
from the fitting. Following convention, we name four quasar 
components A-D, in decreasing order of their brightnesses, 
and we name the lensing galaxy G. The relative positions 
of the four components agree well among the high angular 
resolution images, with a scatter of ~ 0'.'03. The maximum 
separation between images is l'.'76. For the NIR images, the 
lensing galaxy is modelled well by the Sersic profile (n — 4) 
with scale radius Re ~ 0'.'7 — 0'.'9, ellipticity e ~ 0.57, and 
position angle (East of North) 9s ~ 24° . The lensing galaxy 
has somewhat different best-fit scale radius and ellipticity 
in the UH88 V-band image. Re ~ l'.'42 and e ~ 0.28. The 
larger Re in V^-band is in fact consist ent with observations 
of elliptical galaxies (e.g.. [Pahrdl 19991 ) 

In Figure [J] we compare the flux ratios of the quasar 
images derived from the five follow-up images. They are 
broadly consistent between different wavelengths, which fur- 
ther supports the lensing interpretation. However, there are 
some noticeable differences in flux ratios, in particular those 
in the UH88 l/-band. While additional systematic errors as- 
sociated with image fitting might account for the differences, 
a possible physical interpretation is differential dust extinc- 
tion which is com monly seen in lensed quasar systems (e.g., 
iFalco et al.|[l999l ): this helps to explain the different V-band 
flux ratios from those in NIR images, because the effect of 
dust is much more pronounced in V^-band than in NIR. The 
l/-band image was taken 8 — 9 months after the observations 
in near-infrared, and thus time variability (intrinsic and/or 
due to microlensing) of the quasar images might also con- 
tribute to the differences. In particular, the fact that im- 
age B, which is significantly fainter in the V-band image, 
is a saddle-point im age (see which is often de magni- 
fied by microlensing (jSchechter fc Wambsgans3l2002l ) makes 
the microlensing interpretation plausible. Additional imag- 
ing observations in the same and different bands will help 
to distinguish these possibilities. 

We note that the case for gravitational lensing is already 
very strong from its characteristic configuration of image 
components, the similarity of their colors, and the presence 
of a bright lensing galaxy. However, the ultimate confirma- 
tion of its lensing nature will require the spectroscopy of all 
image components to show they have similar spectral fea- 
tures. 

3.2 Optical Spectroscopy 

We conducted spectroscopy of the system with the Dual 
Imaging Spectrograph (DIS) at the ARC3.5m telescope on 
2008 June 6. We used a l'.'5 slit and the gratings of B400 
(dispersion of 1.8A) in the blue channel and R300 (2.3A) 
in the red channel. The slit was aligned to observe nearby 
galaxy G3 (see Figure[l| and SDSS J1330-I-1810 simuffane- 
ously. The spectral resolution is f? ~ 500. We obtained the 
spectrum of the Northern-haff of SDSS J1330-I-1810 (around 
image C and D) where the relative contribution of galaxy G 
to the flux is much stronger; however, because of poor see- 
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Figure 4. Flux ratios of quasar images from the follow-up images 
(see also Figure [3] and Table We plot magnitude differences 
between images B-D and the brightest image A, Am = rriji — mx 
{X = B — D). Dotted horizontal lines indicate median values of 
individual ratios, which we adopt for mass modelling. 
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Figure 5. G-band (left) and Na (right) absorption lines of 
the lensing galaxy at z = 0.373 from the spectrum of 
SDSS J1330-I-1810 obtained with the DIS at the ARC3.5m tele- 
scope. The feature at 5750A is an Mgll doublet at 2: = 1.054 (see 
also Figure O. 



ing (~ l'.'2), we could not separate components C, D, and 
G. The exposure time was 2400 sec. The data were reduced 
using standard IRAF tasks. 

This spectrum shows a weak absorption line at ~ 5900A 
and moderate absorption line at ~ 8IOOA (Figure (5]), in ex- 
cellent agreement with those of the G-band 4304A absorp- 
tion and Na 5889/5896A absorption doublet redshifted to 
z — 0.373, the lens redshift inferred from absorption lines in 
the SDSS spectrum (see Figure [2}. Since these absorptions 
are weak, they were not seen in the SDSS spectrum. We 
note that the DIS spectrum exhibits clear emission lines of 
quasars, which should originate both from images C/D and 
from scattered fluxes from images A/B. 

The redshift of galaxy G3 is 2; = 0.311 (Figure[6l), which 




7000 



Figure 6. The spectrum of galaxy G3 (See Figure[T]l taken with 
the DIS at the ARC3.5m telescope. 



is quite close to the redshift of galaxy G2, z = 0.3126. The 
two galaxies have similar spectra, with significant [Oil] emis- 
sion of the same rest-frame equivalent width. The 4000A 
break of G3 is weaker than that of G2, which implies that 
star formation happened more recently in G3. 



4 ENVIRONMENT OF THE LENS GALAXY 

The wide-field image suggests a possible group of galax- 
ies located near the lens system (see Section [2)| . The high 



noted before jFassnacht & Lubin 20021: Momcheva et al. 


2006: WilUams et al.ll20 


36l:l Auger et al.l 20071: ICabanac et al. 


2007(:iShin et al.ll2008l: 1" 


>eu et al.ll2008h and is theoretically 



the redshift of the possible group, z = 0.31, differs from the 
lens redshift zi = 0.373, a foreground group affects the lens 
potential in a similar way as does a group at the same red- 
shift In this section, we study the distribution and prop- 
erties of galaxies in the field from the SDSS. 

First we extract locations and brightnesses of galax- 
ies in the 12' x 12' field centred on the lens system from 
the S DSS data relea se 6 (DR6 : lAdelman-McCarthv et aP 
l2008h . We adopt the IPetrosianI l| 19761 ') magnitudes in what 
follows. We restrict our analysis to galaxies brighter than 
i — 21, where star-galaxy separation is reliable. To study 
the distribution at z ~ 0.31, we adopt pho tometric red- 
shift measurement in the SDSS dat abases (jCsabai et al.l 
l2003l : lAdelman-McCarthv et al.l 12007*). Specifically, we se- 
lect galaxies whose photometric redshifts are consistent with 
z = 0.31 within their quoted errors. Therefore our study here 
is not necessarily restricted to red elliptical galaxies but in- 
cludes blue galaxies as well. Galaxies with large photometric 
redshift errors, Az > 0.2, are excluded from our analysis. 

We show the angular distributions of all galaxies and 
galaxies at 2: ~ 0.31 in Figure [T] There is a clear concentra- 
tion of galaxies to the southwest of the lens. The structure 
is more pronounced after applying a cut by the photometric 
redshift. We note that Gl is the brightest among galaxies 
at 2 0.31 selected in this way. Together with the agree- 
ment of the spectroscopic redshifts of G2 and G3 (Sections 



^ From the discussion of iKeetonI l l2003t) . it is estimated that a 
weak foreground external shear 7 at 2 = 0.31 is equivalent to an 
external shear of O.87 a,t z = 0.373. 
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Figure 7. Upper: The spatial distribution of all the SDSS galax- 
ies brighter than i = 21 in the vicinity of the lens at (0, 0). North 
is up and East is left. The position of the lens system is indicated 
by a cross. Lower: The spatial distribution of galaxies whose pho- 
tometric redshifts are consistent with z = 0.31. 

[5] and [23]) , we conclude that there is a foreground group of 
galaxies at 2 ~ 0.31 with its centre ~ 120" southwest of the 
lens system (corresponding to transverse physical distances 
of 390/i"^kpc at z = 0.31 and 440/i"^kpc at z = 0.373), cen- 
tred around Gl and G2. It is worth noting that a concentra- 
tion of galaxies can also be seen around G3, which suggests 
that there might be a sub-clump of the group around G3. 

There is a possibility that a concentration of galaxies 
at the lens redshift z = 0.373 exists in additional to the 
group at z — 0.31. However, most of galaxies examined here 
are rather faint and have large errors on the photometric 
redshifts, which prevent us to distinguish structures at z = 
0.31 from z — 0.373. Additional imaging and spectroscopic 
follow-up observations are necessary to explore this issue 
further. 



5 MASS MODELLING 

We constrain the lens potential of this system using the 
observed image positions and flux ratios. The relative po- 




[arcsec] 

Figure 8. The best-fit model (SIEx, flux ratios are included as 
constraints) is shown. See Table [2] for best-fit parameters. The 
black curve and squares indicate the critical curve of the best-fit 
model and the predicted image positions. Observed image posi- 
tions are within the symbols. Grey curves and the square are the 
corresponding caustics and source position. 



sitions of the quasar images and galaxies and their errors 
are adopted from Table [T] Flux ratios are estimated from 
the median of the measurements in the five follow-up im- 
ages (see Figure 131 • Note that the flux ratios agree well with 
those of the T^-band image for which the effect of dust ex- 
tinction is small. Considering the scatter, we assume a con- 
servative error on the magnitude difference between quasar 
images, c7(A?ti) = 0.2. We adopt two standard mass mod- 
els that are widely used in modelling strong lens systems; a 
Singular Isothermal Ellipsoid (SIE) and an SIE plus exter- 
nal shear perturbation (SIEx). Since it is of interest whether 
the models can explain the observed flux ratios or not, for 
each model we consider two cases; one includes flux ratios 
as observational constraints and the other only uses relative 
positions as constraints. The optimisation of model param- 
eters is performed using a software package called glafic (M. 
Oguri, in preparation). 

Best-fit and model parameters are summarised in Ta- 
ble [21 We find that all four models considered here produce 
reasonably good fits to the data (x^/doi < 1). The mod- 
els successfully fit the fiux ratios, therefore this lens system 
does not exhibit any sign of flux ratio anomaly. Moreover, 
the fiux ratios predicted by models with no fiux ratio con- 
straints broadly agree with observed flux ratios. The suc- 
cessful mass modeling further supports the lensing nature 
of this system. We show the critical curve and caustics for 
the best-flt model (SIEx with flux ratios as constraints) in 
Figure [S] which indicates that images B and D are saddle- 
point images and images A and C are minima. Given the 
best-fit mass models, we can predict time delays between 
the images. Assuming that the lensing galaxy is located at 
zi = 0.373, for each best-flt model we compute time delays 
between image A and the other images, which are listed in 
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Table 2. Best-fit mass models of SDSS J1330+1810. The column 'Flux' shows whether flux ratios are included as constraints or not. 
The parameter (?Eini and 7 denote the Einstein radius, ellipticity, and external shear, respectively. The position angle of ellipticity 
and shear, and are measured East of North. Xposi Xgap ^'^d xlux' indicate chi-square values from image positions, lens galaxy 
positions, and flux ratios. The total chi-square is the sum of these three. Also shown are time delays between images predicted by the 
best-fit models, adopting the lens redshift z; = 0.373. 



Model 


Flux 


x'/dof 


Xpos 


A, gal 


Xflux 


^Ein 


e 


ee[deg] 


7 




AtAB [day] 


At AC [day] 


At AD [day] 


SIE 


No 


2.10/3 


2.10 


0.01 




0"97 


0.31 


34 






-0.15 


6.04 


-11.68 


SIE 


Yes 


7.28/6 


2.80 


0.01 
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Figure 9. The UH88 V-band image around the lens system. 
North is up and East is left. Shear directions in our best-fit mod- 
els {6-Y = —88° for the model without flux constraints, and —75° 
for the model with flux constraints) are shown by thick solid lines. 
Directions to a nearby bright galaxy (G3) and the centre of the 
group (G1/G2) are indicated by thin solid lines. 



Table |2l Time delays are rather model dependent, partic- 
ularly for At AS and AtAC, for which the maximum frac- 
tional differences between different models are > 70%. This 
is because shorter time delays are more ea sily affected b y 
perturbations in the primary lens potential (|Ogurill2007al ). 

The large ellipticities of 0.31 — 0.56 in the best-fit mod- 
els are broadly consistent with the observed shape of the 
lensing galaxy (see Section l3.ip . We find that the best-fit 
position angle of the galaxy, 9e , is different between models 
with and without external shear: Oe ^ 34° when the external 
shear is not included, whereas 9e ~ 25° for models includ- 
ing external shear. The observed position angle of galaxy G, 
de = 24° ± 2°, agrees very well with the latter, suggesting 
the non-negligible effect of external she ar on the lens poten - 
tial, as in other quadruple lens systems (|Keeton et al.|[l997l ) . 
However, there is no obvious perturber along the shear di- 
rection, 6.y = —88° or -75° (see Figure O. Thus in Figure 
llOl we present likelihood contours in the 9^-0^ plane. Specifi- 
cally, for each {6e, 6j) we perform minimisations and use 
Ax^ to draw contours at the la and 2a confidence levels. 
We find that the shear direction to galaxy G3 is consistent 
with the observations at the la level. The resultant elliptic- 



-60 



0) 
T3 



-90 



120 



150 



1 1 1 1 1 1 1 1 


' 1 ' 

( 


1 1 1 1 


1 1 






X / 




- 03 








-. (GJ./G2 






1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 


1 1 



10 



20 30 
9, [deg] 



40 



Figure 10. Icr and 2a contours (Ax^ = 2.3 and 6.2 for 2 degree 
of freedom) of the position angles of the ellipticity and external 
shear for the SIEx model. Note that the other model parameters 
are marginalised over. Flux ratios are included as constraints. The 
cross indicates the best-flt values (see Table [2l . The vertical dot- 
ted line indicates the observed orientation of the lensing galaxy. 
Horizontal dotted lines are expected values (see also Figure [9]| 
when the external perturbations are caused by the centre of the 
group (G1/G2) or the nearby bright galaxy (G3). 



ity position angle is 20° < 6*6 ^ 24° and is also consistent 
with the observed orientation of the lensing galaxy. On the 
other hand, the shear direction to galaxies Gl and G2 fails 
to fit the data at the ~ 2a level. Therefore we conclude that 
nearby galaxy G3, rather than the group of galaxies ~ 120" 
southwest of the lens, is likely to be the main perturber to 
the lens system. 

Although the shear amplitudes of 0.05 — 0.1 are slightly 
larger than would be produced by a single galaxy ~ 30" 
away from the lens, the existence of the group suggests that 
dark matter sufficient to produce the shear amplitudes can 
be associated with galaxy G3. This idea is supported by an 
apparent concentration of galaxies at 2; ~ 0.31 around G3 
(see Figure [7|. 

Given the best-fit Einstein radius i?Ein = 0"97, we can 
compute magnitudes of th e lensing galaxy expect ed from 



the Faber-Jacks on relation |Faber fc Jacksonlll976l ). adopt- 
ing the model of lRusin et alT i '2003l ). From the lens redshift 
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zi — 0.373, the apparent magnitudes are predicted to be 
V = 20.7, J = 17.9, H = 16.7 and Ks = 16.1, which 
are 0.7 — 1.2 mag fainter than the observed total magni- 
tudes of the lensing galaxy. Using the same model, we find 
that the observed magnitudes correspond to the lens red- 
shift of zi ~ 0.22. The discrepancy is substantially larger 
than the observed la scatter in the Faber- Jackson relation 
(~ 0.5 mag). While it is not clear what causes the differ- 
ence, a possible explanation is an interaction with nearby 
structure which can strip the outer part of the galaxy mass 
and reduce the Einstein radius without ch anging the overall 
luminosity of the galaxy (e.g., iTreu et"ai . 2008). However, 
there is no obvious companion galaxy at z ~ 0.373 around 
galaxy G, and no obvious signature of interaction in the 
galax;y images. Again, extensive spectroscopy of galaxies in 
the field to study structure around the lens redshift might 
help to resolve this issue. 



6 SUMMARY 

We have presented the discovery of a new four-image lensed 
quasar SDSS J1330-I-1810 {zs = 1.393). This source was se- 
lected as a lens candidate by the SQLS due to its extended 
morphology. Our observations in optical and near-infrared 
indicate that it is a typical fold-type quadruple lens, with 
a maximum separation between images of 1'.'76. From the 
spectrum we measured a lens redshift of zi — 0.373. Stan- 
dard simple elliptical mass models fit the data well, includ- 
ing the flux ratios, implying no evidence for substructure. 
The mass modelling suggests an important contribution of 
the external shear, probably from a nearby bright galaxy. 
There is also a foreground group of galaxies whose centre is 
~ 120" from the lens. The lens galaxy is ~ 1 mag brighter 
than predicted by mass modelling. 

Thus far the SQLS has discovered 31 lensed quasars in- 
cludii ig SDSS J13 30+181CE|, of which 5 are quadruple lenses 
(e.g.. Ilnada et al][2003. : .Kavo et al...2007i ). The low fraction 
of quadruple lenses may imply that the faint end slope of 
the quasar optical luminosity function is shallow, although 
it is ma rginally consis tent with standard theoretical expec- 
tations (|Ogurj|2007"bh . The SQLS has completed ~ 2/3 of 
its survey, implying that a few more quadruple lenses will 
be discovered by the completion of the survey. 
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